Contradiction between strong lensing statistics and a feedback solution to the 

cusp/core problem 
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Standard cosmology has many successes on large scales, but faces some fundamental difficulties 
on small, galactic scales. One such difficulty is the cusp/core problem. High resolution observations 
of the rotation curves for dark matter dominated low surface brightness (LSB) galaxies imply that 
galactic dark matter halos have a density profile with a fiat central core, whereas N-body structure 
formation simulations predict a divergent (cuspy) density profile at the center. It has been proposed 
that this problem can be resolved by stellar feedback driving turbulent gas motion that erases 
the initial cusp. However, strong gravitational lensing prefers a cuspy density profile for galactic 
halos. In this paper, we use the most recent high resolution observations of the rotation curves 
of LSB galaxies to fit the core size as a function of halo mass, and compare the resultant lensing 
probability to the observational results for the well defined combined sample of the Cosmic Lens 
All-Sky Survey (CLASS) and Jodrell Bank/Very Large Array Astrometric Survey (JVAS). The 
lensing probabilities based on such density profiles are too low to match the observed lensing in 
CLASS / JVAS. High baryon densities in the galaxies that dominate the lensing statistics can reconcile 
this discrepancy, but only if they steepen the mass profile rather than making it more shallow. The 
result is contradictory demands upon the effects of baryons on the central mass profiles of galaxies. 
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PACS numbers: 

In the standard cosmological model (known as 
ACDM), the universe is dominated by invisible com- 
ponents called dark energy (A) and cold dark matter 
(CDM). The ACDM cosmology is very successful in ex- 
plaining the cosmic microwave background and the for- 
mation of large scale structure. However, there are chal- 
lenges to ACDM on smaller scales Here we focus on 
the cusp/core problem 0, H, 0] and whether proposed 
solutions to this problem can be consistent with the ob- 
served frequency of strong gravitational lensing. 

One possible solution to the cusp/core problem is tur- 
bulence driven by stellar feedback during galaxy forma- 
tion. If this process drives massive clumps of gas through 
the central regions of the first dark matter halos to form 
[B|, [f|, the central cusp may transform into a soft core. 
Once established, phase space arguments imply that the 
core should persist through subsequent mergers @, Q, 
leading to a final halo profile with a finite core radius for 
all galaxies. Such a situation is consistent with essen- 
tially all kinematic observations [l(| • While this may 
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provide a satisfactory explanation for the cusp/core prob- 
lem, we show here that it leads to a contradiction with 
the mass profiles of elliptical galaxies that dominate the 
statistics of strong gravitational lensing. 

Gravitational lensing provides a powerful tool to de- 
tect dark matter. The lensing efficiency is very sensi- 
tive to the slope 7 of the central mass density profile 
(p oc r~ 7 ). It is well established [ll], [H, EH that when 
galaxies are modeled as a singular isothermal sphere (SIS: 
7 = 2) and galaxy clusters are modeled as a Navarro- 
Frenk- White (NFW: 7 = 1) profile (see Table 1), the 
predicted strong lensing probabilities match the results 
from CLASS/ JVAS. A steeper density slope near the cen- 
ter gives a more efficient lensing rate. For example, if we 
model galaxies with an NFW rather than SIS profile, the 
lensing probabilities are too low compared with observa- 
tions at small image separations [l2| . The presence of a 
central flat core (7 ks 21, m galaxies would further limit 
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the lensing efficiency LJ]. For example, a nonsingular 
truncated isothermal sphere (NTIS), which is an ana- 
lytical model [l6| for the postcollapse equilibrium struc- 
ture of virialized objects, has a soft core that matches 
quite well with the mass profiles of dark matter domi- 
nated LSB galaxies deduced from their observed rotation 
curves. The probabilities for lensing by NTIS halos are 
far too low compared to observations [14j, however. 

In order to investigate the effect of a central core on the 
strong lensing efficiency, we use the density profile of the 
halos directly constrained by observed rotation curves. 
These are well fit [ijj by the cored isothermal sphere 
(CIS). The CIS halo has a finite core radius r c within 
which the density is constant (7 = 0). As well as pro- 
viding a good description of the data, the CIS provides 
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FIG. 1: The rotation curve (left) of the LSB galaxy UGC 5750 
(right). Velocity data come from several independent sources 
and methods, including radio synthesis observations of the 
21 cm spin flip transition of atomic hydrogen two inde- 
pendent [lj|, (2j3] optical long slit observations of the n=3^2 
Balmer transition (Ha), and Densepak integrated field Ha 
spectroscopy [2lj . The various halo types are shown as lines 
(as marked) . The parameters of NFW halos are not free, fol- 
lowing 0, [2l[ from ACDM cosmology. The difference between 
this and the data is the cusp/core problem. The core radius 
of the CIS fit is marked by arrows for the cases of zero and 
maximum disk. For clarity, the full CIS halo is only shown 
for the case of zero disk. Attributing mass to the stellar disk 
detracts from the velocity that can be attributed to dark mat- 
ter (albeit not much in the case of LSB galaxies), increasing 
r c as shown and makes the discrepancy with the NFW pre- 
diction of ACDM more serious. Under no circumstances can 
the halos of LSB galaxies be modeled by SIS. 



a reasonable proxy for the unwieldy NTIS profile. 

The best objects for tracing the mass profile of the 
dominant dark matter component are LSB galaxies. In 
other galaxy types, the stellar mass can provide a non- 
negligible contribution to the rotation velocity at ob- 
served radii. This is not the case for LSB galaxies, whose 
diffuse disks remain dark matter dominated [22| down to 
small radii. These objects persistently suggest that dark 
matter halos possess approximately flat cores 23[ that 
are best fit with CIS halos (Figure [T]) . 

We use the most recent results [24j from a sample of 
LSB galaxies for which rotation curves have been derived 
from high-resolution optical velocity fields. For each halo, 
we calculate the mass M by integrating the CIS density 
profile to the radius r2oo ■ This is the radius of a sphere 
within which the average mass density is 200 times the 
critical density of the universe, typically taken b 2j as the 
virial radius. We c omp ute the halo mass for two brack- 
eting assumptions [24j about the mass of the baryonic 
disk: zero disk, in which the mass of stars and gas is 
neglected, and maximum disk, which attributes the most 
mass possible to the stars without exceeding the observed 
rotation. The primary difference between these two cases 
is in the core radius inferred for the halo. As more mass 
is attributed to the stars, less dark matter is necessary 
at small radii. Consequently, r c grows with stellar mass. 

There are correlations between the core radius r e and 
the halo mass M in both zero disk and maximum disk 
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FIG. 2: The correlation between core radius r c and halo mass 
M. Filled points represent the case of zero disk mass and the 
open points represent maximum disk. The lines are the fit to 
the Densepak data [24| only (circles); fits to long slit data [2(|] 
(squares) give indistinguishable resuts. 



cases. Since our aim is to investigate the effect of the 
core radius on strong gravitational lensing efficiency, we 
fit the relation between r c and M (Figure [5]). As a check, 
we repeat the procedure with independent data [20j . The 
results are indistinguishable. 

The gravitational lensing principle tells us that for any 
spherically symmetric density profile (here, a CIS halo), 
multiple images of a source can be produced only if the 
central convergence k c is larger than unity [251 ]. The 
central convergence is a measure of the central surface 
mass density of the lensing halos. It is both mass and 
redshift dependent. For singular density profiles such as 
SIS and NFW, the central value is divergent, so k c > 1 
is always satisfied and multiple images can be produced 
by any mass. For density profiles with a finite soft core, 
however, the condition n c > 1 imposes a minimum mass 
threshold to produce multiple images. For CIS halos 141 ] . 
we have k c cx M 2 / 3 /r c . The larger the core radius, the 
larger the mass needed to ensure K c > 1. While both the 
zero and maximum disk cases give similar r c -M relations, 



the more conservative case is that with the smaller core 
radius for a given mass; other choices would produce less 
lensing. We thus use the formula fit to the zero disk 
case: r c = 2.25(M/10 12 Af©) 1 / 3 kpc. Interestingly, this 
formula is similar to the one derived analytically in the 
NTIS model QUI]. 

The combined JVAS/CLASS survey forms a well- 
defined statistical sample containing 13 multiply im- 
aged lens systems [3, [53] among 8,958 sources. These 
data provide the image separation A9 for each lens sys- 
tem. The observational probability P bs(> A0) for the 
CLASS/ JVAS survey is shown in Figure [3 

When a remote quasar is lensed by a CIS halo, three 
images are produced. The image nearest the lens is very 
weak. It disappears entirely when the source, lens, and 
observer are aligned, and the Einstein ring appears. The 
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image separation A9 is thus the separation between the 
outer two images. By adopting a model for the density 
profile of lensing halos, their comoving number density, 
and the geometry of the ACDM universe, we can predict 
the properties of the strong lens systems. 

In order to compare with the observed lensing proba- 
bilities, we calculate -Pq[s(> the lensing probability 
for quasars at redshift z s lensed by foreground CIS ha- 
los with image separation larger than A8. The redshift 
z s of the sources (quasars) for the CLASS/ JVAS sam- 
ple has an approximately Gaussian distribution [ll|, [26[ 
with a mean of 1.27 and a dispersion of 0.95. The lens- 
ing rate is sensitive to z s , but the effect of the redshift 
distribution is negligible compared to the choice of halo 
profile. We thus use the mean value z s = 1.27 in our 
calculations. For each lens system, the image separation 
depends on the source position. For the CIS model, how- 
ever, the image separation is almost source position inde- 
pendent [Til ] , so we use the diameter of the Einstein ring 
as the image separation for each lens system. Gravita- 
tional lensing magnifies the brightness of sources, so the 
number of lenses will be overrepresented [28j in any ob- 
served sample. The theoretically predicted lensing proba- 
bility should therefore include a magnification bias (MB) 
correction to the observed probability. The MB is calcu- 
lated on the basis of the total magnification of the outer 
two brighter images [2{|. One of the most important 
elements in predicting lensing probability is the comov- 
ing number density of lensing galaxies. We adopt the 
results recently derived [1(3] from the Sloan Digital Sky 
Survey. The background cosmology is taken from the 
five-year Wilkinson Microwave Anisotropy Probe obser- 
vations [3l]]. The final predicted lensing probability for 
CIS is plotted in Figure [H For comparison, the lensing 
probability of the SIS model is shown with the same pa- 
rameters and approximations as CIS. The NFW model 
14j is also shown. 

The predicted lensing probability for the CIS modeled 
galaxies is about four orders of magnitude lower than the 
observations of CLASS/ JVAS at all image separations, 
and two orders of magnitude lower than the NFW model 
at small image separations. We have used a spherical 
model, as it is kown that the ellipticity does not signifi- 
cantly affect the total lensing efficiency pjjj]. Indeed, the 
combination of all our approximations together can shift 
the result by no more than one order of magnitude, as 
can be seen from the close match of our approximate SIS 
model to the data. So it is safe to conclude that halo 
models like CIS and NTIS with the soft central cores de- 
rived from kinematic observations can not account for 
the statistics of strong gravitational lensing. 

It is not difficult to understand the low lensing prob- 
ability of the CIS model. Recall that the central con- 
vergence depends on the mass M and the redshift Zl 
of the lensing halos. With the fitting formula r c cx 
Af 1 / 3 , we have k c (M,z l ) ~ M 1 / 3 F(zl), where F(z L ) = 



SI{zl) 1/& D l D ls /D s , with Q(z L ) = O ro (l + z L f + n A> 
Dl, Ds and £>ls are the angular diameter distances from 
the observer to the lens, to the source and from the lens 
to the source, respectively. For quasars at zs — 1.27, 
F(zl) has a maximum value of 0.24 for Zi, in the interval 
[0,25]. The condition k c > 1 for strong lensing implies 
M > 3 x 10 13 Af Q . Since M ~ 1O 13 M corresponds to the 
most massive galaxies in the present universe, the galax- 
ies with lower mass provide no contribution to the total 
lensing probability. Furthermore, the contributions of all 
galaxies to the total lensing probabilities are governed by 
the comoving number density n(M), which has a high- 
mass exponential cutoff n{M) ~ exp(-M /3 / 3 ), with 
[3 = 2.67 in our calculations. Consequently, galaxies with 
mass lower than ~ 10 13 Mq make no contribution, and 
high-mass galaxies meet the exponential cutoff. Some 
previous work also used the CIS model for early- 
type galaxies to calculate the strong lensing probabilities, 
and obtained reasonable results. They adopted a typical 
value of the core radius of r c ~ 0.1 kpc, much smaller 
than ours (~ 2.25 kpc), so hardly different from SIS. 

The only difference between CIS and SIS is that CIS 
has a finite core radius. While the SIS model matches the 
lensing observations quite well, the low lensing probabil- 
ities of the CIS model is in serious contradiction to ob- 
servations of strong gravitational lensing. Similarly, the 
NFW/SIS model contradicts rotation curve data. The 
proposed remedy of the cusp/core problem via feed- 
back driven turbulence fixes this problem at the expense 
of creating another. 

Most lensing galaxies are giant elliptical galaxies with 
substantial stellar masses, while we base the CIS model 
on observations of dark matter dominated LSB galaxies. 
These are very different galaxy types. Lensing is not sen- 
sitive to whether the mass doing the lensing is baryonic 
or dark, so the contradiction might be avoided if the total 
mass distribution of ellipticals — stars plus dark matter 
- can be modeled as SIS spheres. The challenge then 
becomes a self-consistent understanding of the formation 
of all galaxy types. 

In the context of the ACDM structure formation 
paradigm, the initial condition for galaxy formation is 
the NFW halo. Baryonic gas dissipates and settles to 
the center of the gravitational potential defined by the 
dark matter to form the visible galaxy. As the gas col- 
lapses, the potential must adjust to the rearrangement of 
mass. This process, commonly referred to as adiabatic 
contraction: 34], has the effect of steepening [3^] the mass 
profile (increasing 7). Since the NFW halo is not ad- 
equate to explain lensing on its own 14, [36| (Figure [3]), 



this process seems necessary to produce elliptical galaxies 
that behave as SIS spheres. Indeed, any transformation 
other than 7 = 1 — > 2 is disallowed. 

In LSB galaxies, we need the opposite process: some- 
thing that drives 7 from 1 — > 0. This is what turbulence 
is proposed^, @] to do. The hypothesized turbulence is 
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FIG. 3: The lensing probability with image separation larger 
than A8. Our prediction for the CIS model based on the ob- 
served r c -M relation (Figure[2]) is shown as the solid line. This 
fails to explain the observed lensing frequency (heavy line) by 
four orders of magnitude. In contrast, our approximate SIS 
model (dotted line) provides a reasonable match to the data. 
A pure NFW model (dashed line) gives intermediate results. 



pealing idea makes predictions that can immediately be 
tested. Elliptical galaxies need to steepen their mass pro- 
file (7 = 1 — ► 2), implying quiescent formation in which 
adiabatic contraction is the dominant effect. Meanwhile, 
LSB galaxies need to experience vigorous bursts of early 
star formation that drive turbulence to erase the initial 
cusp (7 = 1 — > 0). 

It seems quite reasonable to suppose that these very 
different galaxy types formed in different ways. Unfortu- 
nately, all available evidence suggests that these systems 
have behaved in the opposite manner form what we need. 
The stars in ellipticals are predominantly old seem- 
ing to have formed in vigorous early bursts of star forma- 
tion, perhaps associated with frequent merging. They are 
the result of violent, not quiescent galaxy formation. LSB 
galaxies, on the other hand, show no evidence [H, [3{| of 
ever having experience intense star formation. They are 
the poster children for quiescent galaxy formation. This 
is precisely the opposite of what is needed. 

We conclude that the apparent contradiction between 
rotation curves and strong lensing statistics pointed out 
here is genuine. It is difficult to simultaneously reconcile 
the soft cored halos favored by many kinematic obser- 
vations with the singular mass profiles favored by strong 
lensing. In both fundamental tenet of the ACDM 

structure formation paradigm, the NFW halo, is inade- 
quate to explain the observations. Substantial rearrange- 
ment of the initial NFW mass profile is required. Ideas 
hypothesized to solve one problem tend to make the other 
one worse. 
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driven by feedback from early star formation in the first 
halos. If this process is universal and efficient, as pro- 
posed, then we may only solve the cusp/core problem at 
the expense of introducing a new problem with lensing. 
The baryons must first collapse to the center of the halo 
before they can drive feedback there. So only one pro- 
cess can dominate: either adiabatic contraction, which 
increases 7, or feedback, which might reduce 7. If feed- 
back succeeds in establishing a soft core, it should persist 
through subsequent mergers @> HI • It is difficult to see 
how an elliptical galaxy with an SIS mass profile could 
be constructed in this scenario. 

If only some halos experience vigorous star formation 
and feedback while others form their galaxies through a 
gentler adiabatic contraction, we might achieve a consis- 
tent explanation with difference initial conditions lead- 
ing to two very different types of galaxies. This ap- 
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